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INTRODUCTION 


Of 0.6 ,t a MaoK ooefticia’nt 

Laboratory (NPL) 36 transonic wind tunnel (ret. 11 ““‘lonal Physios 

at the British NPL by combinina a "no;iW •• )• e aj.rfoj.i tested was designed 

rear ioadin,. ihe aLrofird:L;.ared LeTpHsro Hh^rr^" Oistribution with 
ot drag coefficient »ith Mach ntnbet. (See ret. 1 I ''>ti»tfon 

creep,** which is defined as a <?irmif' ^ chief prob.lem was **drag 

rise Mach number, data oret ^ a drag occurring below the dLg- 

numbers (2.65 x 10^ to 3.65 x Ld^" ^ obtained at low Reynolds 

in order to determine the aerodynlmic charLteS^L^of 

ues of Reynolds number as well is to determ^n^^r J airfoil at higher val- 

on this type of airfoil, a moLl o^tfe ^^^^510 Lrfo'f condensation 

high Reynolds numbers in the langley 0.3-Meter Transon^ tested at both low and 

TCT). ^ r ^ Meter Transonic Cryogenic Tunnel (0.3-m 


SYMBOLS 

Values are given in SI Units. Measurements 
Customary Units. 

Cp pressure coefficient, ^ 3° 

q 

GO 

^p, sonic pressure coefficient for m = 1 

i 


and calculations were made in U. S. 


model chord, 152.4 


mm 


section chord-force coefficient. 


section 


^ s ^(f) 

profile-drag coefficient, / c' d^— ^ 

wake 


po.'nt drag coefficient (ref. 2) 


section lift coefficient, c^ cos a - c^ sin a 
section pitching-moment coefficient about quarter-chord point, 
-^C^(f-0.3S) ->(h) * ^ C_,(|) d(f) 

section normal-force coefficient, > ^ C d^-j 
vertical height in wake profile, mm 


M 


Mach number 


M local Mach number 

free-stream Mach number 

local static pressure, atm (1 atm = 101.325 kPa) 
total pressure, atm 
free-stream static pressure, atm 
free-stream dynamic pressure, atm 

Reynolds number based on model chord and free-stream conditions 
total tenperature, K 
K airfoil abscissa, mm 

2 airfoil ordinate, mm 

a angle of attack, deg 

Abbreviations: 

AOA angle of attack (see fig. 2(a)) 

0.3-m TCT 0.3-Meter Transonic Cryogenic Tunnel 

NPL National Physics Laboratory 
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APPARATUS AND TESTS 
MDdel 

The wind-tunnel model of the NPL 9510 had a chord of 152.4 mir. (6.0 in.), a span 
of 201.42 mm (7.93 in.), and was constructed of HP 9-4-20 alloy steel. All pressure 
tubing was internal to the model. Fifty flush surface static pressure orifices were 
installed on the model in chordwise rows swept at 9.46*. The upper-surface row con- 
tained 32 orifices and the lower-surface row contained 18 orifices. The orifices 
were connected to a high-precision capacitive potentiometer-type pressure trai.sducer 
(ref. 3) through a scanning valve system. The diauneter of each orifice was 0.254 mm 
(0.010 in.) in order to minimize orifice-induced pressure errors. The coordinates of 
the model are presented in table I, and the shape is shown in figure 1. The model 
orifice locations are contained in table II. In general, the measured coordinates 
are within 0.0254 mm (0.001 in.) of the design values, and the surface finish is in 
the range of 0.102 to 0.203 pm (4 x 10 ^ to 8 x io ^ in.). 
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Langley 0.3-Meter Transonic Cryogenic Tunnel 


The investigation was conducted in the Langley 0.3-Meter Transonic Cryogenic 
Tunnel (0.3-m TCT) with its 203.2- by 609.6-mm two-dimensional test section 
installed. This continuous-flow tunnel uses nitrogen gas as the test medium and can 
be operated over wide ranges of both pressure and temperature to obtain an extensive 
range of Reynolds number. For a 152.4-mm (6.0 in.) chord model, the 0.3-m TCT can 
obtain Reynolds numbers from 1.3 x 10^ at a Mach number of 0.35 to greater than 
45 X 10^ at a Mach number of 0.80. Liquid nitrogen is injected into the tunnel cir- 
cuit to control ten^erature in conjunction with independent control of pressure and 
fan speed to obtain this range of test conditions. Details on the cryogenic concept 
and the operation of the 0.3-m TCT are to be found in references 4 to 10. The valid- 
ity of using cryogenic nitrogen as the test medium has been documented in refer- 
ences 8, 9, 11, and 12. Shown in figure 2 are the test section and other significant 
components of the 0.3-m TCT. 


Wake Rake 

The airfoil drag coefficient was determined by using a wake rake described in 
reference 13 dnd shown in figure 2(c). For the present tests, the rake contained 
five active pitot tubes. As viewed in figure 2(a), pitot tube 1 is 12.7 mm to the 
right of the tunnel midspan. Pitot tube 2 is on the tunnel midspan; tube 3 is 
12.7 mm to the left of midspan; tube 4 is 38.1 mm to the left of midspan, and tube 5 
is 50.8 mm to the left of midspan. Ihe tubes had an outside dicuneter of 1.52 mm 
(0.06 in.) and an inside diameter of 1.02 mm (0.04 in.). The three static tubes of 
the wake rake were not used. However, nine static pressures were measured on the 
sidewall opposite the wake rake. The nine static pressure ports are arranged with 
one port midway between the tunnel floor and the ceiling and with four each, spaced 
25.4 mm apart, above and below this midpoint. Both the pitot and static pressure 
measurements were made in a plane located 112 mm (0.736c) downstream of the trailing 
edge of the airfoil at zero angle of attack. 


Tests and Methods 

The test envelope for the 0.3-Meter Transonic Cryogenic Tunnel and the actual 
test points are depicted in figure 3. The upper boundary of the test envelope is 
fixed by a maximum total pressure of 6 atm and total ten^^eratures (96.6 K to 106.6 K) 
which would give free-stream saturation at a given Mach number. The lower boundary 
is fixed by the maximum total terperature of 327 K at a total pressure of 1.2 atm. 

The Mach number for the present two-dimensional test section is limited to approxi- 
mately 0.9. The filled-in symbol is the two-dimensional equivalent of the cruise 
point based on mean aerodynamic chord for a full-scale present-day wide-body trans- 
port. The flagged symbols at a Reynolds number of 30 x 10^ and Mach number of 0.75 
represent a special condensation effects study to be discussed latter. 

The model had 240-mesh carborundum -grit transition strips in place for all 
but one test at Reynolds numbers less than 6 x 10^. These strips extended from 
4 to 6 percent chord on the upper surface and from 6 to 8 percent chord on the lower 
surface to correspond to the locations of the strips in the NPL tests. The grit 
sizing and transition-strip placement was performed in accordance with the methods 
detailed in reference 14. The transition strips were placed in the following manner. 
The model was masked to expose the area intended for treatment. An adhesive was 
sprayed on the model, and t^e grit was then blown on the treated area to a coverage 
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of 5 to 10 percent. Once dried, the mask was removed, and any adhesive buiid-up that 
occurred along the mask edges was scraped off and polished to deter two-dimensionai 
tripping- 

The airfoil surface pressure distributions measured in the tests were numer- 
ically integrated to generate the normal-force coefficient and the quarter-chord 
pitching-moment coefficient. 

These standard measurements were supplemented with wail static pressure measure- 
ments along both the floor and ceiling of the tunnel. Tbese floor and ceiling mea- 
surements are designed to be used in conjunction with a wall interference correction 
method described in references 15, 16, and 17. The coefficient data presented 
herein, however, have not been corrected. Tube displacement effect correction has 
not been included in the airfoil drag coefficients presented herein. 

PRESENTATION OF RESULTS 

The integrated values of normal-force coefficient c^j, quarter-chord pitching- 
moment coefficient and profile-drag coefficient c^ are given herein, ^cause 

the variation in drag coefficient across the tunnel, as determined wit.i the ftve 
individual pitot tubes, is slight until drag rise is reached, only the drag coeffi- 
cient derived from the midspan tube (tube 2) is presented. 

The uncorrected data are presented as plots of angle of attack a, profile-drag 
coefficient c,, and quarter-chord pitching-moment coefficient c^„ versus normal- 
force coefficient c_. An index to the data figures follows: 


Reynolds 

number. 


1.36 X 10® 


Mach 

number, 

M 

0.3511 


Figure 


1.20 300 






Of 




Reynolds 

number, 

R 

Mach 

number, 

M 

Pt' 

atm 

— 

1 

b' 

< 

Figure 

7.03 X 10® 

0.7143 

2.51 

230 

7(e) 

7.02 

.7338 

2.48 



7(f) 

7.02 

.7527 

2.44 



7(g) 

7.04 

.7729 

2.40 



7(h) 

7.04 

.7925 

2.37 



7(i) 

7.06 

.8096 

2.34 



7( j) 

15.02 

.3977 

5.57 

170 

8(d) 

14.94 

.4975 

4.59 



8(b) 

14.89 

.5968 

3.99 



8(c) 

14.88 

.6933 

3.59 



8(d) 

14.94 

.7133 

3.53 



8(e) 

14.98 

.7336 

3.48 



8(f) 

15.25 

.7536 

3.48 



8(g) 

14.48 

.7727 

3.37 



8(h) 

14.94 

.7921 

3.34 



8(i) 

14.90 

.8112 

3.28 



8( j) 

30.03 

.4993 

5.57 

120 

9(a) 

30.13 

.5987 

4.84 



9(b) 

30.20 

.6970 

4.31 



9(c) 

30.09 

.7145 

4.27 



9(d) 

30.04 

.7351 

5.29 

140 

9(e) 

29.62 

.7575 

2.86 

94 

9(f) 

30.07 

.7561 

5.21 

140 

9(g) 

29.88 

.7743 

4.08 

120 

9(h) 

30.03 

.7933 

4.01 

120 

9(i) 

30.04 

.8155 

3.97 

120 

9( j) 

47.63 

.6990 

5.81 

107 

10(a) 

47.72 

.7188 ! 

5.71 



10(b) 

47.60 

.7366 

5.60 



10(c) 

47.55 

.7580 

5.50 



10(d) 

47.47 

.7763 

5.43 



10(e) 


DRAG CHARACTERISTICS 
Free and Fixed Transition 


The results of an investigation of fixing transition is presented in fig- 
ure 11. Ttie dashed line is a theoretical calculation (ref. 18) which assumes fully 
turbulent flow from the 1 -percent-chord location rearward. These results would indi- 
cate that there is considerable laminar flow at a Reynolds number of 3.60 x io^. In 
addition, it appears that natural transition moves well forward for a Reynolds number 
of 7 X 10^ or at most 10 x 10^. This forward location of transition for the higher 
Reynolds numbers eliminates the need for transition fixing at these values of 
Reynolds number. 
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Effects of Reynolds Number 


Drag polars .- Figure 12 is a summary plot for a Mach number of 0.755. For those 
Reynolds numbers where normal- force coefficient c^ is greater than 0.6, section 
drag coefficient c^j is greater than 0.016. This figure also shows the expected 
decrease in c^ with increase in Reynolds number. For c„ = 0 . 18 , the difference 
in c ^3 between the 7.02 x 10 and 47.55 x 10® conditions is 0.0023. This value is 
comparable to the difference in computed skin frictions for a flat plate at these two 


qrag rise .- Cross-plotting the data of figures 4 to 10 at constant normal-force 
coefficient c^^ and Reynolds number R allows an investigation of the effects of 
Reynolds n^ber on drag-rise Mach number. Examples of this cross-plotting are pre- 
sented in figure 13. If the definition that drag-rise Mach number is the Mach number 


for which 

and below 
cross*-plot 


=0.1 iS applied, the drag rise occurs below M = 0.65 at c^ = 0.35 

M = 0.55 at c^ = 0.55. There were insufficient data at c = 0.55 to 
the data for 47.6 x io ; however, two test points are directly plotted. 


PRESSURE DISTRIBUTIONS 
Conparison of Data 

A conparison of the present data with the NPL data of reference 1 is presented 
in figure 14. The drag values from NPL are slightly greater for c below 0.6 
for the lower Mach numbers (fig. 14(a)) and below 0.5 for the highe? Mach numbers 
(fig. 14(b)). Above c„ of 0.6 and 0.5 for the low and high Mach numbers, resnec- 
present data have higher drag; however, there is variation in the drag 
coefficients measured across the tunnel under these high-lift test conditions. The 
quarter-chord pitching-moment coefficient has very similar characteristics for the 
present and reference data. 


4 -v. . drag-coefficient and angle-of-actack curves (fig. 14), it is apparent 

that the two tunnels (the British NPL 36 and the Langley 0.3-m TCT) have quite dif- 
ferent interaction between the tunnel walls and the model. Ibis difference in wall 
interference is expected since the two tunnels have different geometries. The NPL 36 
funnel has a width of 355.6 mm, a height of 762 mm, and an open area ratio of 0.033. 
The 0.3-m TCT, on the other hand, has a width of 203.2 mm, a height of 609.6 mm, and 
an open area ratio of 0.05. This difference in interaction produces a difference in 
a of about 0.5» for c„ in the range of 0.35 to 0.55. This 0.5» difference is 
confirmed when pressure distributions are conpared, as in figure 15, for the same 
values of c^ and M. in figure 15(a), the normal-force coefficient is 0.44, Mach 
number is 0.5, and Reynolds number is 2.6 x 10®. The reference data have an angle of 
attack of 2.0O, and the present data have an angle of attack of 2.51®. Figure 15(b) 
gives a comparison for normal-force coefficient of 0.42, Mach number of 0.7, and 
Reynolds number of 3.3 x 10 . The present data have an angle of attack of 2.00®, and 
the reference data have an angle of attack of 1.5®. in both cases, only the refer- 
ence data are faired; the agreement in the pressure-distribution results is con- 
sxdered to be very good. 


Comparison of Theory and Experiment 


In figure 16, experimental and theoretical (ref. 17) pressure distributions 
are shown for Mach numbers in the midrange (0.7 to 0.78) and in the upper range 
(0.78 to 0.82). Test points in these two ranges were chosen for comparison because 
they fail very near the 3u x 10^ Reynolds number curve of figure 13(b). The test 
points at a Mach number of 0.736 (fig. 16(a)) have a strong shock well forward on the 
airfoil. The theory of reference 17 has been used to produce the theoretical pres- 
sure distribution by specifying Mach number and normal-force coefficient. The agree- 
ment is rather good considering the fact that shock-wave entropy changes are not 
included in the theoretical conputation. The upper-surface shock strength is under- 
estimated; thus, the region immediately downstream of the shock is also incorrect. 

The test point at a Mach number of C.792 (fig. 16(b)) has a slightly weaker shock 
farther back on the airfoil. In this case, the theory underestimated the strength 
and misplaced the location of the upper-surface shock and overpredicted the extent of 
the supersonic region for the lower surface. These two test points (M = 0.'7i6 
and M == 0.792) were then corrected for blockage by the xerogram of reference 16 to 
obtain the free-air Mach number. These corrected Mach numbers and the measured 
normal-force coefficients were then specified to conpute the pressure distributions. 
The results of these computations are shown in figure 17. 'fhe Mach number has been 
corrected to 0.726 from 0.736, and the angle of attack has been computed as 1.88® 
instead of 1.78® for the midrange case (geometric a being 2.47®). These are small 
corrections which make minimal changes in the conputed pressure distribution. How- 
ever, they are sufficient to inprove the agreement for the lower surface and to 
slightly improve the predicted value of the strength of the upper-surface shock. In 
the upper-range case, the Mach number has been corrected to 0.777 from 0.792, and 
the angle of attack has been conputed as 1.58® instead of 1.40® (geometric a 
being 2.23®). These corrections are slightly larger and make significant changes in 
the conputed pressure distribution. The lower-surface computed distribution is now 
in excellent agreement with the experimental data. The conputed shock strength and 
locatiov are greatly improved for the upper surface. Note that values of both C 
and c^ have been corrected to take into account the corrected free-stream Mach ^ 
number. The theoretical predictions corrected for Mach number and lift coefficient 
give slightly better agreement with the experimental data than the uncorrected 
predictions. 


CONDENSATION 

There is considerable benefit to be made in terms of increased Reynolds number 
capability and reduced operating cost if a cryogenic tunnel is operated as cold as 
possible. (See ref. 11.) This can be illustrated by referring to figure 18. The 
saturation curve labeled M = 0 is the combinations of total pressure p^ and 
total temperature T^ for which nitrogen will normally condense if held at static 
conditions. The curve labeled M_ = 0.7568 defines the conditions for which satu- 
ration would occur in the test section without a model present and with a free-stream 
Mach number of 0.7568. The curve labeled M =1.5 defines the conditions for 
which saturation would occur when a local Mach number of 1.5 is reached on the model. 
The dashed line represents all of the possible combinations of total pressure and 
total temperature that yield a Reynolds number of 30 x 10^ based on a 152.4-mm chord 
for this free-stream Mach number of 0.7568. The test point at 94 K and 2.86 atm 
(square symbol) is very close to free-stream saturation, and condensation may occur 
on the model in sufficient amounts to affect the aerodynamic characteristics. On the 
other hand, the test point at 140 K and 5.21 atm (circle symbol) is chosen such that 
a local Mach number of 2.18 on the model would be needed for condensation to even be 
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possible. This latter test point requires 2.22 times as much liquid nitrogen to 
maintain as does the former test point at 94 K and 2.06 atm. For these two sets of 
test conditions, the maximum local Mach number on the model was 1.5, and there are no 

apparent condensation effects for 94 K as indicated by the aerodynamic data presented 
in figures 19 and 20. 

The significance of operating as cold as possible at the same Reynolds and Mach 
numbers without altering the aerodynamic data is an immediate reduction in operating 
cost. If it is possible to operate as close to the free-stream saturation curve 
(temperature wise) as the 94 K condition - but at higher pressure - without altering 
the aerodynamic data, the gain is in Reynolds number. For example, 101.5 K and 

5.21 atm give a Reynolds number of 50 x 10^ at a Mach number of 0.7568 (triangle 
symbol in fig. 18). 


CONCLUSIONS 

A set of aerodynamic data has been presented for an NPL 9510 airfoil. The test 
Reynolds number ranged from 1.34 x 10^ to 48.23 x io^. Significant conclusions 

foil nw - 


1. The drag creep detected by the British National Physics Laboratory (NPL) 

during previous tests at low Reynolds numbers is also present at high 
Reynolds numbers. 

2. The section drag coefficient continued to decrease with Reynolds number to 

the highest Reynolds number covered during these tests. 

3. The model was tested very close to free-stream saturation without altering 

the aerodynamic data by condensation effects. 

4- The present tests produced the same pressure distribution as the previous NPL 
tests when normal-force coefficient, Mach number, and Reynolds number were 
matched. The geometric angle of attack between the previous and the 
present tests is different by about 0.5®. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 19, 1983 
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TABLE 1.- AIRFOIL COORDINATES 
Upper Surface 


x/c 

z/c 

.0000 

-.00002 

.0004 

.00392 

.0016 

.00760 

.0025 

.00942 

.0050 

.01277 

.0100 

.01670 

.0150 

.01922 

.0200 

.02100 

.0250 

.02238 

.0300 

,02357 

.0400 

.02553 

.0500 

.02722 

.0600 

.02871 

.0700 

.03007 

.0800 

.03129 

.0900 

.03241 

.1000 

.03343 

.1200 

.03530 

.1400 

.03694 

.1600 

.03839 

.1800 

.03970 

.2000 

.04090 

.2200 

.04205 

.2400 

.04315 

.2600 

.04417 

.2800 

.04507 

.3000 

.04592 

.3200 

.04669 

.3400 

.04741 

.3600 

.04806 

.3800 

.04867 


x/c 

z/c 

.4000 

.04917 

.4200 

.04959 

.4400 

.04990 

.4600 

.05013 

.4800 

.05027 

.5000 

.05031 

.5200 

.05028 

.5400 

.05016 

.5600 

.04995 

.5800 

.04970 

.6000 

.04930 

.6200 

.04879 

.6400 

.04818 

.6600 

.04741 

.6800 

.04655 

.7000 

.04543 

.7200 

.04415 

.7400 1 

.04263 

.7600 

.04093 

.7800 

.03902 

.8000 

.03688 

.8200 

.03459 

.8400 

.03205 

.8600 

.02933 

.8800 

.02642 

.9000 

,02337 

.9200 

.02009 

.9400 

.01655 

.9600 

.01278 

.9800 

.00888 

1 .0000 

.00490 




TABLE I AIRFOIL COORDINATES - Concluded 


Lower Surfoce 



z/c 

.6800 

-.01705 

.7199 

-.01064 

.7600 

-.00503 

.8000 

-.00047 

.84-00 

.00299 

.8800 

.00475 

.9200 

.00485 

.9600 

.00322 

1 .0000 

.00002 


x/c 

z/c 

.0000 

-.00002 

.0004 

-.00397 

.0016 

-.00795 

.0025 

-.00976 

.0050 

-.01353 

.0100 

-.01846 

.0150 

-.02201 

.0200 

-.02498 

.0250 

-.02753 

.0300 

-.02984 

.0400 

-.03388 

.0500 

-.03738 

.0600 

-.04043 

.0700 

-.04310 

.0800 

-.04547 

.0900 

-.04753 

.1000 

-.04940 

.1200 

-.05265 

.1600 

-.05759 

.2000 

-.06114 

.2400 

-.06340 

.2800 

-.06438 

.3200 

-.06408 

.3607 

-.06247 

.4000 

-.05962 

.4400 

-.05553 

.4800 

-.05036 

.5200 

-.04428 

.5600 

-.03762 

.6000 

-.03067 

.6400 

-.02373 
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TABLE n.“ ORIFICE LOCATIONS 


Upper Surface 


Lower Surface 


x/c 

z/c 

.00000 

.00000 

.00997 

-.01833 

.01768 

-.02362 

.05263 

-.03820 

.10228 

-.04976 

.15266 

-.05678 

.20203 

-.06125 

.27703 

-.06435 

.37565 

-.06150 

.45068 

-.05426 

.52567 1 

-.04337 

.60065 ! 

-.03065 

.67547 

-.01 770 

.71730 

-.01103 

.85073 

.00358 

.90095 

.00503 
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Figure 2.- Test section and significant corrponents of Langley 0.3-Meter 

Transonic Cryogenic Tunnel. 






Model turntable 



(b) Schematic drawing showing major cortponents 
Figure 2.- Continued. 




Concluded 




Figure 3.- The 0.3-m TCT test envelope and test points. (Flagged symbols represent 

special condensation effects study.) 
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Figure 5.- Continued, 










Figure 6.- Section characteristics for two different values of Reynolds number 
number, total pressure, and total tenperature with free transition. 






Figure ?•- Continued 












characteristics at Reynolds numbers 
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